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^y. , Abstract 

pj ■ We construct the formal solution of the Cauchy problem for the 

dispersionless Kadomtsev - Petviashvili equation as application of the 
Inverse Scattering Transform for the vector field corresponding to a 
Newtonian particle in a time-dependent potential. This is in full anal- 
ogy with the Cauchy problem for the Kadomtsev - Petviashvili equa- 
^r^ \ tion, associated with the Inverse Scattering Transform of the time 

CN ' dependent Schrodinger operator for a quantum particle in a time- 

Z^ ■ dependent potential. 

O 

\0 • 1. Dispersionless (or quasi-classical) limits of integrable partial differential 

equations (PDEs) arise in various problems of Mathematical Physics and are 

intensively studied in the recent literature (see, f.i., [B 121 El 13 El)- In par- 

^ \ ticular, a quasi-classical dressing has been developed |11 for the prototypical 

example of the dispersionless Kadomtsev - Petviashvili (dKP) (or Khokhlov- 

>. Zabolotskaya) equation: 

rS 

c3 : Utx + Uyy + {uu^)^ = Q, u = u{x,y,t) eM, x,y,teR. (1) 

In this paper we construct the formal solution of the Cauchy problem on 
the plane for the following system of PDEs in 2-|-l dimensions: 

Uxt + Uyy = -{uux)x - VxUxy + VyUxx, u,v eR. x,y,t e M, . . 

^xt ~r ^yy ^ ^xx ^x ^xy ~r ^y ^xx 

and for its f = reduction, the dKP equation (|T]), as application of the 
recently developed Inverse Scattering Transform (1ST) for vector fields [0]. 



O 



Indeed the system (0) arises as the compatibihty condition of the Lax pair 

Li^ = 0, L2^ = 0, (3) 

implying [Li,L2] = 0, where Li, L2 are the following vector fields: 
Li = dy + {p + v^)dx - u^dp, 



2 , ^. , . ^. (4) 



L2 = dt + {p^ + pv^ + u- Vy)d^ + {-pu^ + Uy)dp. 

Setting ti = in (j3)), one obtains the Lax pair of the dKP equation, which 
was derived in [Sj taking the quasi-classical limit of the well-known Lax pair 
of the KP equation [ZUHj. 

We remark that, in the dKP reduction v = Q, the two vector fields are 
Hamiltonian and the Lax pair Q takes the form 

i^y + Pi^x - U^i^p = i'y + {Hi, '0}(p,a;) = 0, ,^. 

A + (P^ + U)i^x + {-PUX + Uy)'4'p = A + {H2, ^}{p,x) =0, 

in terms of the two Hamiltonians [S] 

2 3 

Hi = — + u, H2 = — +pu- d~^Uy, (6) 

where {■, ■}{p,x) is the standard Poisson bracket with respect to the canonical 
variables {p, x): 

{/, 9}{p,x) = fp9x - fx9p, (7) 

leading to the Hamiltonian form of dKP: Hn — H2y + {H2, i^i}(p^a;) = 0. 

Since the Lax pair Q of the dKP-like system (j21) is made of vector fields, 
Hamiltonian in the dKP reduction (^, the eigenfunctions satisfy the follow- 
ing basic properties. 

1) The space of eigenfunctions is a ring. If /i, /2 are two solutions of the 
Lax pair Q, then an arbitrary different iable function F(/i, /2) of them is a 
solution of 0. 

2) In the dKP reduction v = 0, the space of eigenfunctions is also a Lie 
algebra, whose Lie bracket is the natural Poisson bracket ^. If /i, /2 are 
two solutions of the Lax pair ©, then their Poisson bracket {fi, f2}(p,x) is 
also a solution of (0). 

2. Now we consider the Cauchy problem for the dKP system (j2)) and for the 
dKP equation (^ within the class of rapidly decreasing real potentials m, v. 

M,v^O, {x^ + y'^)^oo, ue^, (x,?y)GM^ t > 0, (8) 



interpreting t as time and the other two variables x, y as space variables. To 
solve such a Cauchy problem by the 1ST method [^, we construct the 1ST 
for the operator Li, within the class of rapidly decreasing real potentials, 
interpreting the operator L2 as the time operator. 

The localization (jHJ of the potentials u, v implies that, if / is a solution 
of Lif = 0, then 

f{x,y,p) ^ f±{^,p), y^±oo, ,g. 

^■= x-py; 

i.e., asymptotically, / is an arbitrary function oi ^ = x — py and p. 

A central role in the theory is played by the two real Jost eigenf unctions 
V^i,2(3^; y,p), the solutions of I/iV?i,2 = uniquely defined by the asymptotics 

^i{x,y,p) ^ ^, (P2ix,y,p) ^ p, y ^ -00. (10) 

In this paper we often use the compact vector notation: / = (/i, f2)^- Then: 

The Jost eigenfunction is the solution of the linear integral equations = 
^ + G{—Vx0x + Ux0p), for the Green's function G{x, y,p) = 6{y)S{x — py). 
The y = +00 limit of defines the natural scattering vector a for Li: 

hm ip{x, y,p) = S{0 = e + <?(0. (12) 

J/— »+oo 

The direct problem is the transformation from the real potentials u, v, 
functions of the two real variables {x,y), to the two real scattering data 
o"i,cr2, the components of the scattering vector a, functions of the two real 
variables (^,p). Therefore the mapping is consistent. The impact of the dKP 
reduction t) = on these and other data will be shown below. 

A crucial role in the 1ST theory for the vector field Li is also played by 
the analytic eigenfunctions ip±{x,y,p), the solutions of Lii{j± = satisfying 
the integral equations 

'if±ix, y,p) = J^2 dx'dy'G±{x -x',y- y',p)[-Vx'ix', y')ilj±^,{x' , y',p)+ 
Ux' (x', y')ip±pix', y', p)] + f, 

(13) 



where G± are the analytic Green's functions 

G±ix,y,p) = ±—- , . X 1 - (14) 

27n[x — [p ± ie)y\ 

The analyticity properties of G±{x, y,p) in the complex p - plane imply that 
ip+{x,y,p) and ijj_{x,y,p) are analytic, respectively, in the upper and lower 
halves of the p - plane, with the following asymptotics, for large p: 



Mx,y,p) = ^ + lUix,y) + O i^^j , \p\ »1, 

-yu{x,y) -v{x,y) \ (1^) 



It is important to remark that the analytic Green's functions ()14|) exhibit 
the following asymptotics for y -^ ±oo: 

entailing that the y = +cxd asymptotics of'ip+ and ip- are analytic respectively 
in the lower and upper halves of the complex plane C,, while the y = — oo 
asymptotics of 11)+ and ip- are analytic respectively in the upper and lower 
halves of the complex plane ^ (similar features have been observed first in 
PH] and later in 0). 

The Jost eigenfunctions Lpi^2 form a basis; thus any solution / of Li/ = 
is a function of ip. The analytic eigenfunctions 'il)± possess the representations: 

^± = /C±((^) = (^+ x±(<^), (17) 

defining the spectral data x±- 

Since the y — > — oo limit of (fTTj) reads: 



hm 4-?"*=X±(6, (18) 



y- 



the above analyticity properties of the LHS of (fTH|l in the complex ^ - plane 
imply that x+{^ and X-H) ^'^^ analytic respectively in the upper and lower 
halves of the complex plane ^, decaying at ^ ~ oo like 0(^~^). Therefore 
their Fourier transforms x+{^) ^i-nd X-{^) have support respectively on the 
positive and negative oji semi-axes. 



The spectral vectors x± can be constructed from the scattering vector a 
through the following linear integral equations 



(19) 



x+{^) + ^(^i) y^{^) + /k2 dr] x+iv)Qiv, ^) j = 0, 

X^{uJ) + 9{-uJi) (^a{uj) + /jg2 df]X-{f])Q{ff,uj)j = 0, 



involving the Fourier transforms a and x± of f? and x±'- 

~a{^) = [ d^a{Oe-'^<, ^±(cD) = / d^x±iOe"^-^^ (20) 

and the kernel: 

Q{ff,co) = [ JL^e^iv^-^Xyv-^^iO _ 1]. (21) 

To prove this result, one first evaluates (fTTj) at y = +cxd, obtaining 



( hm 4 - i) = m + X±(e + '?(0). 



(22) 



Applying the integral operator f^^ d^e~'^'^'^- for cui > and cji < respec- 
tively to equations (|12I)+ and (J22I)-, using the above analyticity properties 
and the Fourier representations of x± and a, one obtains equations (jT^. 
The reality of the potentials: u,v & M. implies that, for p G M, (^ = (^, 

'?/'+ = "0-; consequently: a = a, x+ = X-- 

3. An inverse problem can be constructed from equations (jTTj). Subtract- 
ing ^ from equations (fT7|) - and (fT7j) +. applying respectively the analyticity 
projectors P+ and P_: 

P =±— f ^^P^_ f23) 

^ 27rzyi,p'-(p±ze)- ^ ^ 

and adding up the resulting equations, one obtains the following nonlinear 
integral equation for the Jost eigenfunction 0: 

^i /r 7^h^)X+{'f{x, y,p')) = C 



Once (f is reconstructed from (j24j) . the analytic eigenf unctions follow from 
(J17p . and u,v from equation (J15p . This inversion procedure was first intro- 
duced in [TT] and also used in j6j. 

4. As u,v evolve in time according to ^, the t-dependence of the spectral 
data S and /C±, defined in (IT^ and (fT7|). is described by the equations: 



Ei(e,p,t) = t(S2(e-p't,p,o)r + si(e-p%p,o) 



(25) 



where Si and S2 are the two components of the vector S, identifiable with 
each of the spectral vectors S and /C±. To prove it, we first observe that 

(j)i{x,y,t,p) = ipi{x,y,t,p)-tipl{x,y,t,p), .^g) 

(I)2ix,y,t,p) = ip2ix,y,t,p) 

are a basis of common Jost eigenfunctions of Li and L2. The y = +00 limit 
of equation 1/202 = yields S2t + P^S2^ = 0, while the y = +00 limit of 
equation 1/201 = yields {dt +p^(9^)(iSi — fiS|) = 0, whose solutions are (j^ 
for S. Analogously, 

are a basis of common analytic eigenfunctions of Li and L2; therefore 

7r±i = IC±i{<pi, 02), vr±2 = /C±2(01, 02), (28) 

for some functions IC±i^2 depending on x,y,t,p only through 0. Comparing 
at t = these equations with equations (fT7|) . one expresses /C±i,2 in terms of 
/C±i 25 obtaining equations (j^ for /C±^ 2- 

We observe the unusual resonant character of the explicit t-dependence 
(J2S1) of the spectral data, if compared to the more elementary one, obtained 
in ;^, for the heavenly equation [T^ . 

5. In the Hamiltonian dKP reduction v = 0, the transformations ^ — > S{^), 
C, -^ ^±{0 ^^re constrained to be canonical: 

{5i,52}^-={/C±i,/C±2}5-=l. (29) 

To prove it, we observe that the Poisson bracket of the eigenfunctions (pi 
and (/92 is also an eigenf unction: Lp^ = {'Pi,'P2}{x,p), Lip^ = 0. Using 



the asymptotics (fTUI). one infers that tp^ ^ 1, at y ^ —oo; therefore, by 
uniqueness, Lp^ = 1. Evaluating now the Poisson bracket (f^ at y = +00 and 
using (fT^ . one obtains the constraint (j^ for S. We also observe that the 
eigenfunctions {ijj+^,ijj+2}(x,p) and {ilj^^,ilj-2}ix,p) are analytic in the upper 
and lower p plane and go to 1 at \p\ -^ 00. Since 1 is also an eigenfunction, 
by uniqueness they are identically 1: {i/j^^, 11)^2} {x,p) = 1- Therefore, from 
the equations: 

{^±l,^±2}{x,p) = {l^±l, !<^±2}{V1,'P2){VI, ^2}{x,p) = 1, (30) 

consequence of (fT7|) . one infers the constraints (j^ for /C±. 

6. It is well-known (see, f.i., jl3j) that linear first order PDEs like (jHj).(PJ) 
are intimately related to systems of ordinary differential equations describing 
their characteristic curves. The Hamiltonian dynamical systems associated 
with the vector fields L12 of dKP are: 



dx 



^=p'^ + U,= {H2,x}(^p^^), 

J = -PU^ + Uy= {H2,p}(p^x), 



(32) 



Therefore the dKP equation characterizes the class of time - dependent po- 
tentials for wich the Newtonian flow ^Sl\) commutes with a flow with cubic, 
in the momentum p, Hamiltonian. 

There is also a deep connection between the above 1ST and the y-scattering 
theory for the commuting flows (jHTj) and (jH^ . Let 0(x, y, t,p) be the basis of 
common eigenfunctions of Li and L2 defined in fl26|) : then, solving the system 
uj = (j){x,y,t,p) with respect to x and p (assuming local invertibility), one 
obtains the following common solution of (|^ and (jH^ : 

^ r/ , \ f ^ \ ^/ , ^\ ( ^2V + ^"^i + ^i\ 

u = (p{x,y,t,p) ^ ] =r{y,t;uj) r^ / ],y 00. 



P J \ ^2 

(33) 
The y = +00 limit of the solution f{y, t; uj): 



defines the scattering vector A{u) = fi{uj) — Co oi (|M|) and (jH^ . which is 
connected to the 1ST data S by inverting the system u = S{x — py — p'^t,p, 0) 
with respect to x and p: 

a = S(.-py-pH,p.O) «(;) = ( "=(-)!' +^i(^)' + "■(-) ) . (35) 

The transformation Co — > Q{Co) is clearly canonical: {^i, f22}(^i,(^2) = 1. 

Since the dynamical system (j31|) describes the motion of a Newtonian par- 
ticle in the plane subjected to a generic time - dependent potential u{x, y), as 
a byproduct of the 1ST of this paper one can reconstruct, from the scattering 
vector A{Co) of the dynamical system (jSU, the time dependent potential u. 

Remark 1. The are two other ways to do the inverse problem. The first one is 
the linear version of the nonlinear problem (j24j) . obtained exponentiating the 
Jost and analytic eigenfunctions used so far. Consider the following scalar 
functions: 

$(x,2/,p; a) = e^°-^(^*P), '^±{x,y,p; a) = e*"-^^±(^'2^'P), a G M^ (36) 

Due to the ring property of the space of eigenfunctions, also $(x,y,p;a) 
and "^ ±{x , y , p; a) are eigenfunctions; ^{x,y,p;d) is characterized by the 
asymptotics $ -^ exp{ia ■ ^), y -^ — oo, while '^±{x,y,p;a) are analytic 
respectively in the upper and lower halves of the p plane, with asymptotics: 
^^ = exp{id-i)[l+p-^a- U{x,y) + 0{p-'^)]. 

Exponentiating the representations (fT7|) . one obtains the expansions of 
the analytic eigenfunctions \E'± in terms of the Jost eigenfunction $: 

^±ix,y,p;d) = ^{x,y,p;d) + J^2d(3K±{d,(3)^{x,y,p;f]), 

K^idJ) ^ /,, 4,e'("-^>«le-^-x±«) - 1]. ^'^^ 

Multiplying the equations ^7\ + and ((HZI)- by exp{—id ■ ^), subtracting 1, 
applying respectively P_ and P+, and adding the resulting equations, one 
obtains the following linear integral equation for $: 

(38) 



in which we have omitted, for siinphcity, the parametric dependence on (x, y). 
Once $ is reconstructed from (jHS|) and, via (jHZI), ^± are also known, the 
potentials are reconstructed in the usual way from the asymptotics of \l/±. 

The third version of the inverse problem is a more traditional (nonlinear) 
Riemann-Hilbert (RH) problem. Solving the algebraic system p7|) - with 
respect to Lp: tp = L{ip_) (assuming local invertibility) and replacing this 
expression in the algebraic system (fT7j)-|-. one obtains the representation of 
the analytic eigenfunction ^+ in terms of the analytic eigenfunction -f^.: 

V^+ = 7?(V^_) =^5*- + W-), peM, (39) 

which defines a vector nonlinear RH problem on the realp axis. The RH data 
TZ are therefore constructed from the data /C by algebraic manipulation. 
Viceversa, given the RH data R, one constructs the solutions %jj± of the 
nonlinear RH problem (|39|) and, via the asymptotics (|15p. the potentials. 

— * 

As for the other spectral data, one can show that the t-dependence of 7^ 
is described by (j^ and the dKP constraint reads {TZi,7l2}(^,p) = 1, while 

the reality constraint takes the form: 1Z{1Z{^, ^)i^) = ^, V(^, for p G M. 
Remark 2. Dressing schemes can be formulated from the three different 
inverse problems presented in this paper in a straightforward way. 
Remark 3. The 1ST constructed in this paper allows one to solve the Cauchy 
problem for the whole hierarchy of PDEs arising from the commutativity 
equation [I/i,I/2 ] = 0, where the coefficients of the vector field Lg are 
polynomials in p of arbitrary degree n &N. 

Remark 4- There are deep similarities between the Cauchy problem for dKP 
and the Cauchy problem for the heavenly equation, recently solved in [B|, 
since they are both based on the 1ST for Hamiltonian vector fields (the 
dKP equation is actually a geometric reduction of the heavenly equation 
|H]). There is, however, an important difference between these two cases. The 
vector fields of the dKP equation contain partial derivatives with respect to 
the spectral parameter p, unlike the case of the heavenly equation jB] . 
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